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We determine the curvature of the phase transition Hne in the /i — T plane through an analysis 
of various observables, including the Polyakov loop, the quark number susceptibilities and the 
susceptibility of the chiral condensate. The second derivative of these quantities with respect to 
jU was calculated. The measurements were carried out on A^t- =4,6,8 and 10 lattices generated 
with a Symanzik improved gauge and stout- hnk improved 2+1 flavour staggered fermion action 
using physical quark masses. 
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1. Introduction 

Quantum chromodynamics (QCD) is the theory of strong interactions. While at low temper- 
atures the dominant degrees of freedom are hadrons, due to one of the most important properties 
of QCD, asymptotic freedom, at high temperatures it describes a different phase of matter called 
quark-gluon plasma (QGP). The phase transition between the hadronic phase of matter and QGP 
can be investigated by lattice simulations. The transition at zero chemical potential - which repre- 
sents the case of equal number of quarks and antiquarks - is of huge importance, since it is relevant 
for the early Universe and also for high energy collisions. The region of the phase diagram corre- 
sponding to small and moderate pL is also interesting, since the cooling down of the QGP after a 
heavy ion collision occurs in this area. 

Simulations at nonvanishing chemical potential are burdened by a significant problem: the 
fermion determinant here becomes complex, and as a result makes importance sampling impossi- 
ble. A possible solution for this issue is to expand observables into a Taylor-series in /i, where the 
coefficients can be calculated at vanishing chemical potential. Due to the symmetry of the partition 
function describing the system, the first term in these expansions always vanishes. The second term 
on the other hand is related to the curvature of the transition line. 

The 2-1-1 flavour QCD transition was found to be an analytic crossover ^ (instead of a first- 
order phase transition), which usually results in different transition temperatures for different ob- 
servables [Q]. Our aim is to determine the curvature k: of the phase transition line for different 
observables. These include the Polyakov loop, the strange quark number susceptibility, the chiral 
condensate and the chiral susceptibility. 

2. Possible scenarios 
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Figure 1: Two possible scenarios about the QCD phase diagram around zero jX. Transition temperatures 
determined by different observables may converge to (left side, compare to first lattice indications on Nj — 
4 [^) or diverge from (right side, first lattice indications on Nt —4- [|||) each other as the chemical potential 
is increased. This tendency can also have an effect on the existence (and if it exists, the position) of the QCD 
critical endpoint. 
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As an effect of the crossover character of the transition at zero jx, the critical temperature 
determined from the chiral susceptibility {Tc{x^/\it) ~ 151 MeV) and that from the strange quark 
number susceptibility {Tc{Xs) ~ 175 MeV) are quite different [Q]. This situation can change at finite 
chemical potential according to two possible scenarios. The strengthening of the transition could 
cause these two values to approach each other, while if the transition remains a weak crossover, 
the two transition lines can diverge (see figure |l|). The two cases are separated by the sign of the 
difference k:(Xs) — ^{Xw) curvatures for the two observables. A positive sign (left side) 

would indicate that the crossover region shrinks, which might suggest that a critical endpoint exists 
on the phase diagram. On the other hand a negative sign would cause the crossover region to 
expand, which may result in the abscence of the critical point. 

Supposing that there exists a critical endpoint around baryonic chemical potential /ig ^ 360 
MeV, the difference of the curvatures should be around Ak f» 0.02 (see definition below). Here 
and also in the following under /i we mean the /ig baryonic chemical potential. There are several 
works in the literature in which the value of the curvature was computed [|, g 0, g |]. In these 
Nt =4 works the curvature was found to be in the interval K = 0.003 . . . 0.01. On the other hand, 
the Nt = 6 results of [||] indicate that the curvature might become less pronounced towards the 
continuum limit. 



3. The curvature 



Let us parametrize the transition line in the vicinity of the vertical jJ. = axis as Tc{lJ,) 
Tc (l — K" • n^/T^) . This implies that the curvature can be written as 



d(M^ 



(3.1) 



We illustrate our first procedure to measure the derivative in question using the quantity Xs/T'^ 
(see definition in section |^. Since < XsjT'^ < 1 stands independently of \i, we can define rc(/i) 
via the relation Xs/T'^{Tc,}Jl) = 0.5. In the following we use the short notation ^ = Xs/l"^- 

The total derivative of the observable i^(r,/i^) may be written as di^ = {d&/dT) - AT + 
(d / d {pL^)) -dpL^. So, since along the Tc{}Jl) line = by definition, for the derivative in the 
definition of the curvature this reads as 




(3.2) 



This way it is enough to measure these two derivatives at an arbitrary temperature. For more 
complicated quantities this procedure gives only an approximate value for the curvature. In these 
cases one should carry out measurements at various temperatures; then (/x) can be defined as 
the inflection point (for e.g. Xs/T'^) or the maximum value (for e.g. X^iw^ of the observable. This 
second procedure - which is clearly more CPU-demanding - is an ongoing project. Nevertheless, 
to calculate the derivative d& /d{}JL^) we need to measure new operators. These are detailed in the 
next section. 
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4. Technique 

Let us consider the partition function in its usual form 

^ = J me-^>^^"\detMff^^ (4.1) 

and denote the derivative with respect to /!„ j by '. The derivatives of 3f are easily calculated to be 
(log = {riu.d) and (log = {Xu,d), where the light quark number density «„ ^ and the light 
quark number susceptibility Xu,d are the following combinations: 



nu4 = ^Tr {M-'m') 



Using these definitions the second derivative of any (possibly /i„^ -dependent) observable can be 
determined as 



2 - \- Au,d/ - \^/\Zu,d) + {2^'nu.d + (4.2) 



d^ld 

For observables that do not depend explicitly on iJ.u,d (like the Polyakov loop L or Xs) the third 
term in ( |4.2| ) vanishes. For -dependent observables the derivatives ff' and were calculated 
numerically, using a purely imaginary chemical potential. 



5. Simulation setup 

We used a Symanzik improved gauge and stout-link improved staggered fermionic lattice ac- 
tion in order to reduce taste violation. The configurations were generated with an exact RHMC 
algorithm. We determined the line of constant physics (LCP) using physical masses for the light 
quarks m„ as well as for the strange quark m^. The LCP was fixed by setting the ratio tiik/ fx and 
niK/mji to their physical values. We used four different lattice spacings Nt = 4,6,8, 10 and aspect 
ratios Ns/Nj of 4 and 3. The scale was fixed by /k and its unambiguity checked by calculating niK, 
fjc and tq. For measuring the operators necessary for the above mentioned derivative the random 
noise estimator method was used. The number of random vectors was set to 80, in order for the 
error coming from the random estimator method and that from the finiteness of the statistics to be 
of the same extent. The details of the simulation setup can be found in [Q] or |[To|]. 



6. Observables 



In order to determine the derivative dTc/d{jJL^) expressed in (3.2), one needs to select some 
observable For this role we used the following quantities: 
The Polyakov loop is defined as 

Y Nt-i 

L=-3£Tr n ^4(^,0 (6.1) 
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In order to be able to extract continuum limit results, an appropriate renormalization is necessary. 
For L this means L,- = Lexp(V(ro)/2r), where V{r) is the T = static potential 
The strange susceptibility 

needs no renormalization, since it is connected to a conserved current. 

The chiral condensate can also be expressed as a derivative of the partition function: 

r^iogiT 

W=T7^ (6-3) 

V dm 

It can be renormalized by subtracting the additive divergences, and then multiplying by the quark 
mass, so that the multiplicative factors also cancel (here the factor is used to get a dimensionless 
combination): 

V/i/A,-= (i//vA-i//v^(r = 0))-m-^ (6.4) 
Finally, the definition of the chiral susceptibility is 



_ r^^iogjT 
Cij/i// — 

and the renormalization is the following 



~ T7 5 2 (6-5) 



7 1 

Xwr = {Xw-XwwiT = 0))-m^-j^ (6.6) 



7. Results 



As already mentioned, in our approach we select an arbitrary temperature to calculate the cur- 
vature of the transition line. In order to increase statistics, we can also average measurements at 
different temperatures. This procedure is illustrated on figure |^ for the case of the strange suscep- 
tibility and the chiral condensate. 

Having obtained the renormalized quantities for Nj =4,6,8 and 10, we are in a position to 
carry out the continuum extrapolation. The results shown on figure ^ are only preliminary; note 
that the extrapolated values have quite large statistical errors. 



8. Summary 

While the curvature values for 1//!//^ are smaller than those for Xs/T^ every lattice spacing, 
they become consistent in the a — > limit. Our aim in the first place is to determine the difference 
K"(i//i//r) — k{Xs/T^), so increasing the statisctics in order to clarify the situation is very important. 
Nevertheless we report the continuum results: K{\jfYr) = 0.0034(40) and k{Xs/T^) = 0.0013(51), 
and also for the Polyakov loop and the chiral susceptibility: K{Lr) = —0.0095(93) and K{x^i\iir) = 
-0.0018(34). 

As the next step we are looking forward to perform the above described procedure to measure 
dff/d{ll^) for the whole temperature interval around T^. This way it becomes possible to determine 
the critical temperature for the given observables (by locating their inflection point or maximum 
point) at nonvanishing chemical potential. Using this technique and much larger statistics we plan 
to give a final answer for the curvature of the QCD transition line. 
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Figure 2: The curvature determined using the normalized strange susceptibility (left side) and the renor- 
malized chiral condensate (right side), respectively. Independent measurements at different temperatures (j3 
values) can be averaged to decrease statistical errors. Results from A^j = 8 lattices are shown. 




Figure 3: The curvature of the transition Une as a function of the lattice spacing squared, for the strange 
susceptibility (left side) and the chiral condensate (right side). As characteristic to the applied action, the 
A^r = 4 results are out of the scaling regime, so the linear extrapolation is carried out using only the finer 
lattices. As it can be seen, larger statistics is needed in order to extract a reliable continuum result. 
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